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q ■ Abstract 

We study cosmic-rays in decaying dark matter scenario, assuming that the 
dark matter is the lightest superparticle and it decays through a i?-parity 
violating operator. We calculate the fluxes of cosmic-rays from the decay of 
the dark matter and those from the standard astrophysical phenomena in the 
same propagation model using the GALPROP package. We reevaluate the 
preferred parameters characterizing standard astrophysical cosmic-ray sources 
with taking account of the effects of dark matter decay. We show that, if 
energetic leptons are produced by the decay of the dark matter, the fluxes 
of cosmic-ray positron and electron can be in good agreements with both 
PAMELA and Fermi-LAT data in wide parameter region. It is also discussed 
that, in the case where sizable number of hadrons are also produced by the 
decay of the dark matter, the mass of the dark matter is constrained to be 
less than 200—300 GeV in order to avoid the overproduction of anti-proton. 
We also show that the cosmic 7-ray flux can be consistent with the results of 
Fermi-LAT observation if the mass of the dark matter is smaller than ~ 4 TeV. 



1 Introduction 



An anomaly has been discovered in the recent observation of cosmic-ray positron. 
The PAMELA satellite has observed positron excess in the flux around the energy 
range of 10—100 GeV, which cannot be explained by standard astrophysical phenom- 
ena in our Galaxy [lj. A possible solution is to consider non-standard phenomena 
in astrophysics such as nearby pulsars [2]. On the other hand, the PAMELA exper- 
iment has motivated us to settle this problem in the viewpoint of particle physics; 
the anomaly can be explained by the decay or the annihilation of dark matter of the 
universe. Indeed, various scenarios with decaying EH — 1.211 and annihilating [2~2~] — [2Uj 
dark matter have been discussed in much literature]^ 1 

A well-motivated scenarios is decaying dark matter scenario in the framework of 
supersymmetric model. Supersymmetric model provides natural candidate for dark 
matter, i.e., the lightest superparticle (LSP). In addition, supersymmetric model is 
one of the most attractive models beyond the standard model; it gives a solution 
to hierarchy problem and enables the three gauge coupling constants to unify at 
~ 10 16 GeV. A very long lifetime of the LSP, which is required to solve the PAMELA 
anomaly, can be realized with a very weak i?-parity violation (RPV). In the scenario 
with unstable LSP dark matter, not only positron but also electron, 7-ray, and/or 
anti-proton give significant contribution to cosmic-ray fluxes. Thus, in order to solve 
the PAMELA anomaly in this framework, it is necessary to check if the fluxes of 
all the cosmic-ray particles are in agreement with observation. For such a purpose, 
one should calculate the fluxes of cosmic-rays from standard astrophysical source, 
i.e., supernova remnant (SNR), as well as those from the decay of dark matter. In 
many of the previous studies, however, fluxes of cosmic-rays from different sources 
are evaluated using different propagation models and/or different SNR injection 
parameters. For a complete study of the decaying dark matter scenario, the fluxes 
of all the cosmic-ray species should be calculated in a consistent way. 

In this paper, we consider the case that dark matter is unstable (and long-lived) 
and it becomes the significant source of high energy cosmic-rays. We pay particular 
attention to the case where the LSP is dark matter, and study the cosmic-ray fluxes 
in light of PAMELA and other cosmic-ray observations. We use the GALPROP 

^For the decaying dark matter scenario, it was pointed out that unstable dark matter seems 
to be suggested by the HEAT results [30] and that the signal of decaying dark matter could be 
observed by PAMELA as a positron excess before the PAMELA result showed up [HJ 132] • 
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package [33J so that the cosmic-ray fluxes from different sources are evaluated with 



the same propagation modelscj We reevaluate some of the parameters character- 
izing the SNR injection spectra taking account of the effects of the decay of dark 
matter. We then discuss constraints on the decaying scenario from observations. 
We will show that, by properly choosing model parameters, the PAMELA anomaly 
can be well explained in the decaying dark matter scenario without conflicting other 
constraints. 

The organization of this article is as follows. In the next section, we explain how 
we calculate cosmic-rays from the standard astrophysical phenomena and decaying 
dark matter. Here important parameters in our analysis are described. Numerical 
results are then shown in section [3j Section 0] is devoted to conclusion. 



2 Cosmic-Rays from SNR and Dark Matter 

Spectra of cosmic-ray particles which we observe depend on what the sources of 
cosmic-rays are and how the particles from the sources propagate. In this section, 
we first consider cosmic-rays from SNR which is supposed to be standard source of 
high energy cosmic-rays. Cosmic-rays from decaying dark matter are then discussed. 
In our numerical calculation, all the cosmic-ray fluxes are obtained by using the single 
numerical code, GALPROP, in order to perform a consistent analysis. We explain 
relevant parameters characterizing energetic cosmic-rays from astrophysical sources 
and dark matter decay. 



2.1 Cosmic-rays in standard Galactic model 

Cosmic-rays mainly consist of nuclei, electron, positron, anti-proton, and 7-ray. 
Among those, electron and nuclei, e.g., proton, helium, carbon, oxygen, and iron, 
are considered to be from the remnants of supernovae and pouring to the earth after 
they have drifted by interaction with interstellar matters and magnetic field in our 
Galaxy. These are called primary cosmic-rays. On the other hand, cosmic-rays are 
also generated secondarily in a consequence of interaction processes; we call them 
secondary cosmic-rays. (See Table [TJ) In the collisions of cosmic-ray with interstellar 
gas, such as pp-collision, anti-protons and other nuclei such as lithium, beryllium, 
boron, and sub-Fe (scandium, titanium, and vanadium), are produced. In the col- 

# 2 For the former study of decaying scenario by using GALPROP, see Refs [4], |34 ) — [40 ] . 
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Table 1: Propagation and production of cosmic-ray e^, p, p, and 7 from SNR. Here e prir 
(Pprim) and e^ c (p sec , psec) are primary and secondary e~ (p) and (p, p), respectively. 



lision processes, pion and kaon are also produced; secondary positron and electron 
are then emitted in the cascade decay of such mesons. 

The primary spectra of electron and proton from SNR are assumed to obey power 
law, since they are assumed to be produced through Fermi-acceleration mechanism. 
The source terms of these cosmic-ray particles are then parametrized as 

QsNR = A--E , GeV' (2-1) 
QsNR = A»PGeV (2-2) 

Here, A e - and y e (A p and j p ) are normalization and power index of electron (proton) 
injected from SNR, while -EceV (PGev) is GeV-normalized energy (momentum) of 
electron (proton). Following the treatment adopted in "Conventional model" of 
cosmic-ray [IT], W2\. the power indices are assumed to take different values for high 
and low energy region as in Table [2j If SNR is the dominant source of high energy 
cosmic-rays, the shapes of those spectra are sensitive to the index parameters 7 e and 
7 P . In particular, primary electron generally dominates the total (e + + e~) flux, and 
its spectrum in the energy region over 10 GeV is sensitive to the value of 7 e . 

If we adopt the Conventional model, predicted cosmic-ray fluxes such as boron to 
carbon (B/C) ratio, proton, helium, and anti-proton fluxes are in good agreements 
with observations. The values of parameters adopted in the model are summarized in 



Table I21F 3 I Within this model, however, cosmic-ray e + from SNR cannot explain the 
positron fraction data reported by PAMELA, which is now known as the "PAMELA 
anomaly". In addition, the (e + +e~) spectrum observed by Fermi-LAT [13] is slightly 



* 3 One can consider diffusive reacceleration or diffusive convection; we adopt the former one. 
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Table 2: Important parameters of cosmic-ray used to calculate fluxes from SNR. The pri- 
mary electron (proton) is normalized to fit observations, and j e (7 P ) is chosen to reproduce 
the observed spectrum. Those values are adopted in Conventional model. 



harder than the prediction of the model; however, as we will discuss later, such a 
discrepancy may be solved by adopting an appropriate value of 7 e . 

In the study of high energy cosmic-ray, it is also important to consider cosmic 7- 
ray and radiation fluxes. Inside the Galaxy, high energy 7-ray is necessarily produced 
by several processes (see Table [1]): the inverse- Compton (IC) scattering of e ± with 
interstellar radiation field (ISRF), the bremsstrahlung of with interstellar gas, the 
synchrotron radiation from e ± under magnetic field, and the decay of pion which is 
produced in the hadronic reaction of cosmic nuclei in interstellar gas. The spectra 
of 7-ray and synchrotron radiation depend on the injection spectra of SNR electron 
and proton. According to Ref. jH], the Conventional model gives 7-ray flux well 
below the observed flux and radiation flux consistent with observation. It is also 
shown that the consistency is hold if 7 s < 2.0 for E < 10 GeV is satisfied [4"T] P 4 
Electron in this energy range is, however, considered to be subject to relatively large 
effect of the solar modulation. 

Since we are interested in the scenario where cosmic-rays are produced from the 
decay of dark matter, we do not necessarily restrict ourselves to the Conventional 
model. We should rather reevaluate some of the parameters characterizing the SNR 
injection spectra with taking account of the effects of the decaying dark matter. In 
this paper, we take A e - and 7 e to be free parameters. (The other parameters are 
taken as the same as in Conventional model unless otherwise mentioned.) 

Hereafter, the cosmic-rays originating in SNR are called as background (BG): 



$ e iG = *e". +$e- c , (2-3) 

BG prim scc 

= $ e+ , (2-4) 



# 4 In order to be consistent with local e flux in 1 GeV < E < 30 GeV, the power index of 
injected electron j e ~ 2.5 — 2.7 for E > 0(1 GeV) is favored. 
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where $ - and $ e ± are fluxes of primary e and secondary e ± , respectively Be- 

prim sec 

sides electron and positron, hadrons are also produced which result in background 
proton and anti-proton. We also calculate BG fluxes of these particles: 

^PBG = ^Pprim + ^Psec) (2-5) 

$pbg = ^W, (2.6) 

where $p prim is primary p flux, while $ Psoc ($p soc ) is secondary p (p) flux. In the 
calculation of cosmic-ray fluxes of p and p, we use force-field model [H] with the 
use of the solar modulation potential of <fi = 550 MV|^j In addition, for later 
convenience, we introduce normalization parameter (donated as a e ), which is de- 
fined by $ e - (E) = a e &°} (E), with being the reference flux normalized as 

P r ^ m ^prim ^prim 

$ (0) (E = 100 GeV) = 1.3 x lO" 11 GeV^nrtec^str- 1 . 

^prim 

For the estimation of the BG 7-ray flux, we also use the GALPROP package. 
There may exist other astrophysical contributions which are not taken into account 
in the GALPROP package; we do not consider such^contributions in our analysis 
because they are expected to have large uncertainties 



#6 



2.2 Cosmic-rays from decaying dark matter 

Now we consider cosmic-rays from the decay of dark matter. In the decaying scenario, 
it is assumed that the lifetime of dark matter is much longer than the present age of 
the universe. Most of dark matter therefore survives today. Even though the decay 
of dark matter is suppressed by the long lifetime, it can be a significant source of 
cosmic-ray if standard-model particles are produced through the decay. Cosmic-rays 
from dark matter depend on the spectra of particles injected from the decay and 
how they propagate in the universe. 

Important quantities in the calculation of comic-ray fluxes from unstable dark 
matter are lifetime of dark matter (r DM ), mass of dark matter (m DM ), and spectra 
of emitted particles. With these, source term of cosmic-ray is given by 

Qdm - (2.7) 
tdm^dm dE 



* 5 This is a method proposed in order to take into account the effect of solar modulation. Esti- 
mation of the effect of solar modulation is, however, still uncertain. We thus focus on the energy 
region where the effect is small enough in our numerical analysis. (See the discussion below.) 

#6 Unidentified cosmic 7-ray may have various origins, for example, galaxy clusters [45) . ener- 
getic particles in the shock waves associated with large-scale cosmological structure formation [46 , 
distant gamma-ray burst events, baryon-antibaryon annihilation |47j . 
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where X is cosmic-ray particle {e.g., X = e ± ,p,p, and 7), and dN x /dE is the spec- 
trum of X from the decay of a single dark matter. In addition, p^M is mass density 
of dark matter. In the decaying scenario, the fluxes of cosmic-rays are insensitive to 
the profile of dark matter density except for comic-ray 7 from Galactic center. Since 
we will not study such 7-ray, we adopt the isothermal profile for Galactic halo; 



where p® ~ 0.43 GeV/cm 3 is the local halo density, r core ~ 2.8 kpc is the core radius, 
r Q ~ 8.5 kpc is the distance between the Galactic center and the solar system, and r 
is the distance from the Galactic center. We have checked that our numerical results 
are almost unchanged even if we use other dark matter profiles such as the NFW 
profile [4"8] . 

In the decaying dark matter scenario, cosmic 7-ray from extra- Galactic region 
may have significant flux because the 7-ray is produced by the IC process induced 
by the energetic e 1 * 1 from dark matter decay [H]. It was shown that the flux could 
be comparable to the observed 7-ray flux. (For the discussion of constraint from 
isotropic 7-ray observation, see also [501 EI]-) We calculate the extra-Galactic 7-ray 
by following formula given in [49]. In the calculation of cosmic 7-ray from extra- 
Galactic region, we use averaged dark matter density aa #7 



Notice that the 7-ray from the extra-Galactic region is isotropic. 7-ray is also pro- 
duced at the central region of our Galaxy. However, the flux from the Galactic center 
strongly depend on the dark matter profile [50]. Thus, in order to derive conservative 
constraint, we only consider 7-ray from high Galactic latitude. 

In supersymmetric model, there are many candidates for dark matter which can 
decays under RPV. One of the possibilities motivated by the PAMELA anomaly is 
unstable dark matter which decays into final state which consists of only energetic 
leptons. Then, charged leptons can make prominent rise in positron fraction to 
explain the anomaly without conflicting the observation of cosmic-ray anti-proton. 
In such a case, 7-ray is inevitably produced by IC process as secondary cosmic-ray. 
(Those production mechanism are summarized in Table [31 including the following 

* 7 Because we only consider the 7-ray flux averaged over large solid angle in the decaying dark 
matter scenario, we do not have to worry about the effects of sub-halo. These effects may be 
important in the annihilating dark matter scenario as shown in Refs. |52j — [58]. 





(extra-Galaxy) 

Pdm 



1.2 x 10~ 6 GeV/cm 3 . 



(2.9) 
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Final state 



Primary Secondary 



Leptons — > e — > 71c 

Leptons + hadrons — > e^, p, p, 7 — > 71c 

Table 3: Production process of cosmic-rays in the decaying dark matter scenario we 
consider. Here, 71c means 7-ray produced in inverse-Compton scattering process. 



case to be described below.) On the other hand, one can also consider a case where 
significant amount of hadrons (as well as leptons) are contained in the final state of 
the decay of dark matter. In such a case, anti-proton as well as 7-ray may provide 
stringent constraints on the scenario. Here, notice that scenario with similar final 
state results in similar cosmic-ray fluxes. We thus focus on several cases which 
give typical results; gravitino and sneutrino dark matter scenarios with bi-linear 
and/or tri-linear RPV. In each case, we calculate dN x /dE, including contributions 
of particles form cascading decay after hadronization, with the use of the PYTHIA 
package [59] . 

Let us first consider the case where gravitino is the LSP. In the present case, 
because of the smallness of RPV, the production mechanism of gravitino is unaffected 
by the RPV interaction; gravitino can be produced by scattering processes of thermal 
particles [60] . decay of the LSP in minimal supersymmetric standard model sector 



(MSSM-LSP) after freeze-out [60, 61J or a non-thermal process [62jcf| (Sneutrino, 
which we will discuss below, can be also produced by non-thermal process |63j). 
The spectra of final-state particles, on the other hand, strongly depend on how the 
i?-parity is violated. Because we are interested in the PAMELA anomaly, we only 
consider RPV operators by which the LSP decays into final states with energetic 
lepton. Then, one possibility is to introduce the following soft breaking bi-linear 
RPV operators: 

£ bi _RP V = B l l l H u + m\. H liH* d + h.c, (2.10) 

where Lj is left-handed slepton doublet (with % = 1 — 3 being the generation index) , 
while H u and Hd are up- and down-type Higgs boson doublets, respectively. (Here, 
Li and Hd are defined in the bases in which the bi-linear RPV terms in superpotential 
are rotated away.) With the operators given in Eq. ( I2.10p . gravitino decays into two- 



* 8 However, the production via the decay of the MSSM-LSP is suppressed if the lifetime of the 
MSSM-LSP becomes shorter than the cosmic time of its freeze-out. 
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body final states: — l v h ^ v ii Wh> and hui. (Here and thereafter, gravitino is 
denoted as ip^.) Decay width of each process is given in Ref. [3T]; the decay process 
is dominated by the mode ip^ — > WU if the gravitino is heavier than iy-boson. The 
lifetime is related to the bi-linear RPV parameter as [3 lj , 



where m 3 / 2 is gravitino mass. In addition, n 2 = K h where 

Bi sin (3 + m\ cos (3 
m = ^ , (2.12) 

with tan (3 = (H®) / (H%) (here (H°) and (H%) are the vacuum expectation values of 
up- and down- type Higgs bosons, respectively), and mo i being the mass of sneutrino 
in z-th generation. We will see that the lifetime of dark matter should be O(10 26 sec) 
in order to explain the PAMELA anomaly. Such a lifetime is realized when k ~ 

o(io- 9 - io- 10 ). 

Another important possibility is the decay through tri-linear RPV: 

W-toi-Rpv = -XijkLiLjE^., (2.13) 

where Lj = (0 Li , l Li ) and Ef are left-handed lepton doublet and right-handed lepton 
singlet, respectively, where "hat" is for superfield. With the operator given in Eq. 
(12.131) . gravitino decays as ip^ — > v^-F^y. through diagrams with virtual slepton. 
Here, for simplicity, we assume that the right-handed sleptons are much lighter than 
the left-handed ones. The energy distribution of the final-state leptons is given 

tflV-K/iiifl _ A 2 m 3/2 4(1 - z R ) 



dE lL dE lR 768^71^ (mr /m 3/2 ) 2 - (1 - z R ) ' 



(2.14) 



where Mpi ~ 2.4 x 10 GeV is the reduced Planck scale, mj is the right-handed 
slepton mass, and Zr = 2Ei R /m 3 / 2 . In our analysis, Ay*; is determined so that 
the preferred value of the lifetime of gravitino is obtained. Notice that, taking 
rrif ~ m 3 /2 for simplicity, the lifetime is estimated as 

^.Tx^secx^)- 2 ^)" 3 . (2,5) 

Here, Xijk must satisfy the condition A^ < 10~ 7 in order not to wash out baryon 
asymmetry [64J; this constraint is satisfied in the parameter space we will be inter- 
ested in. 
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Table 4: Scenarios of decaying dark matter considered in this paper. Model parameters 
giving the "best-fit" results to cosmic-ray observations are also shown in each case. 



If the sneutrino is the LSP, monochromatic leptons are produced by its decay 
via the tri-linear RPV superpotential given in Eq. (I2.13p . Indeed, the sneutrino can 
decay as v — > I^Ir- Because the sneutrino is a viable candidate for dark matter 
irrespective of its handedness [65j E3], and also because the fluxes of high energy 
cosmic-rays are sensitive to the spectra of emitted particles from the decay of dark 
matter, we also consider the case of unstable sneutrino LSP. In this case, the lifetime 
is determined as [66] . 

o m26 Asin^V 2 / m f \-i 

t d ~ 2 x 10 sec — — , (2.16) 

V !0" 26 / V200 GeVJ ' v ; 

where 9^ is the sneutrino mixing angle. 

3 Numerical Results 

Now we are at the position to quantitatively discuss the cosmic-ray fluxes. As we 
have already mentioned in the previous section, we assume that the dominant sources 
of the cosmic-rays are supernova remnants and decaying dark matter. In order to 
study the cosmic-ray fluxes from both sources using the same propagation model as 
well as to take account of the production of secondary cosmic-rays, we utilize the 
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GALPROP package for the calculation of all the cosmic-ray fluxes (except for the 
extra-Galactic 7-ray flux). 

The spectral shape of each cosmic-ray particle from dark matter depends on 
properties of the decay of dark matter. We take m^M and r DM as free parameters in 
this paper, and consider three typical cases summarized in Table. HJ the cases with 
(I) the emission of a charged lepton as well as a weak boson resulting in energetic 
jets after the decay (i.e., the gravitino LSP with bi-linear RPV), (II) the emission 
of only non-monochromatic leptons (i.e., the gravitino LSP with tri-linear RPV), 
and (III) the emission of only monochromatic leptons (i.e., the sneutrino LSP with 
tri-linear RPV). It should be noted that there are other candidates for the LSP; the 
lightest neutralino is a popular one. Cosmic-ray fluxes from the decay of the lightest 
neutralino are similar to those obtained in the gravitino LSP case. We thus omit to 
study the case of neutralino LSP. 

For the BG electron flux which originates in SNR, we take 7 e and the normaliza- 
tion of SNR electrons to be free parameters. We then discuss when the decaying dark 
matter scenario gives consistent result with the present cosmic-ray observations. At 
the beginning of this section, we first give our numerical procedure to determine the 
"best-fit" values of the model parameters. Then, we will show the results of our 
numerical analysis. 



3.1 Numerical procedure 

In order to determine the values of the model parameters preferred by the results of 
PAMELA and Fermi-LAT experiments, we define the x^-waiiable as 

2 — 2,2 /o 1 \ 

Xe XPamela ' X Fermi - \ / 

Here, variables Xp am eia an d XFermi are defined as 

N P ( fl( th ) _ R(f 8 A 2 N F ( $( th ) _ . _ $(° bs ) 



.-> \ y J 2 v y e++e ,i e^+e ,i j 

XPamola — / 4 A/? 2 ' Xftermi — / j ' 

where R^+\ (R^+f) is simulated (observed) positron fraction in the i-ih bin, and 
^e++e- i (^i++e- J ^ s s i mma t e d (observed) total (e + +e~) flux. AR e + i and A$ e + +e - ti 
are observational errors in the i-th bin, and A^p and Np are the number of data points 
for those. Since effects of the solar modulation to cosmic-ray e ± is conspicuous in the 
energy range below O(10 GeV), we use five data points of the PAMELA experiment 
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above 15 GeV in the calculation of Xpameia- Then, Np = 5 (whereas Np = 26). In 
addition, we neglect the systematic error from the energy calibration in the Fermi- 
LAT data unless otherwise mentioned; the effect of the systematic error is not so 
important for most of the cases. 

We calculate the xl~ variable defined above as a function of m DM , tdm, 7 e , and 
a e . Because our primary purpose is to find a solution to the PAMELA anomaly, we 
first calculate Xp am eia with mp,u being fixed, and identify 95% C.L. allowed region 
(corresponding to Xpameia < H-07) on the (tdm, 7 e , a e )-space. After that, in the 
parameter region allowed by the PAMELA data, we find out parameter space which 
satisfies x\ < 44.99, which corresponds to the 95% C.L. allowed region (for 31 degrees 
of freedom). We also search the point where the total x\ * s minimized; we call such 
a point as the "best-fit" point. Using the parameters, we also simulate cosmic 7-ray 
and cosmic-ray p and p to check consistency with those latest observations. 



3.2 Gravitino dark matter with bi-linear RPV 

To begin with, we study the case where the gravitino dark matter decays dominantly 
into the first generation lepton through the bi-linear RPV operator, which is the case 
when B\ 3> B2, B%. Using the xt~ variable, we have found the 95% allowed region 
at m 3 /2 ~ 200 GeV and m 3 / 2 > 1 TeV. The best-fit parameters in this case are 
summarized in Table @] as cases (1-1) to (1-3), where the values of Xpameia an d xl 
in each case are also shown. Positron fraction and the total (e + + e~) flux are 
depicted in Fig. [1] (top two panels), where observational data are also given in each 
panel; the PAMELA data [1J on the top left panel, while Fermi-LAT [43], ATIC [67], 
PPB-BETS [68J, and HEAT data [30J on the top right one. 

Our result shows that, if m 3 / 2 ~ 200 GeV, the decaying gravitino LSP may 
explain the PAMELA data with being consistent with the Fermi-LAT observation. 
In the present monochromatic is emitted in the decay, resulting in a sharp 

edge in the positron fraction at E ~ 100 GeV. On the other hand, in the total flux, e 
from dark matter is overwhelmed by BG flux, and hence the edge in the total flux is 
almost hidden by BG. Even though a small edge is visible in the total flux, the value 
of Xe indicates that the model is allowed at 95 % C.L.. (See Table HI) In addition, 
we have also checked that the XFeaaL,^ 1 ^ 8 ^ ^ calculated solely from the total flux is 
smaller than the 95 % C.L. bounded Thus, we conclude that reasonable agreements 



# 9 It is also described in [BS] that the scenario which gives such a sharp edge in the flux is not 
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Figure 1: Cosmic-rays from the decaying gravitino LSP in cases (1-1) to (1-4) in Table HI 
where positron fraction (top left), total (e + + e~) flux (top right), anti-proton flux (middle 
left), anti-proton to proton ratio (middle right), extra-Galactic 7-ray flux (bottom left), 
and the flux from inside and outside of the Galaxy (bottom right) are shown. In each 
panel, results with 7 P = 2.42 (corresponding to case (1-1), (1-2), and (1-3), where gravitino 
mass is taken as m 3 / 2 = 200 GeV, 1 TeV, and 4 TeV, respectively) are shown in magenta 
lines from left to right, whereas blue line is used for those with 7 P = 2.52 (corresponding 
to the case (1-4) where m 3/ / 2 = 200 GeV is chosen). Signal from the decaying dark matter 
and astrophysical BG are depicted by dotted and dashed lines, respectively, while solid 
lines correspond to signal + BG. 
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with both the PAMELA and Fermi-LAT data are realized when 777,3/2 ~ 200 GeV. 
Such a possibility has not been e xplic itly described in previous works after the Fermi- 



LAT (e + + e~) data was released P 10 ! On the other hand, with larger gravitino mass, 
agreement with the PAMELA data becomes good but not with the Fermi-LAT data. 

Concerning cosmic-ray proton and anti-proton, the p flux is shown in the middle 
left panel of Fig. [TJ Here we take the same model parameters as in top two panels. 
Observational data by BESS [70], CAPRICE [7T] and PAMELA [72J experiments are 
also shown in the panel. It can be seen that the signal from the decay of dark matter 
is sm aller than background in most of the energy range, irrespective of the gravitino 
mass P 11 ! However, if we compare the simulated result with the latest observation 
of the ratio between anti-proton and proton (p/p) by the PAMELA experiment, the 
scenario turns out to be strongly constrained. On the middle right panel in Fig. [TJ 
we plot the ratio, where observational data of the PAMELA experiment [72T ITT] is 
also shown. From the figure, it can be seen that the ratio becomes significantly 
larger than the observed one if the gravitino mass is large. Even in the case of 
777.3/2 ~ 200 GeV, the predicted ratio is slightly larger than the observation. 

Here, we should recall that the BG proton flux depends on injection index and 
normalization of SNR nucleon spectrum, whose values are determined by fitting the 
observed flux of proton, B/C ratio, and so on. By taking account of the uncertainties 
of these observations, one can vary the values of injection index and normalization 
these parameters, which may relax the severe constraint. Let us take 7 P = 2.52 
instead of 7 P = 2.42 above 4 GeV with the n orm alization being unchanged from 
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The result with this choice is 



the value adopted in the Conventional model 
summarized in Table H] as the case (1-4), and p flux and p/p ratio are depicted in 
middle two panels of FigfTJ (blue line). We calculated \ 2 variable based on the p/p 
ratio given in the latest PAMELA data [72J. Then, we found Xp/ P = 12.3 (with 
the degree of freedom being 23), and hence the p/p ratio becomes consistent with 



excluded by \ 2 analysis. 

# 10 This possibility in the framework of the decaying dark matter with a few hundred GeV mass 

was pointed out in the scenario where the dark matter decays dominantly into pair |35j . 

* n Some people may think that the result seems to be inconsistent with the previous understanding 

about p from the hadronic decay. In fact, the flux simulated in the transport equation within the 
GALPROP code gives that between MED and MEN models proposed in Ref. [73] . 
# 12 By the explicit simulation of proton and B/C ratio, we have checked that their spectra are 
not affected by this change over 10 GeV. Less than 10 GeV, the flux is significantly affected by 
solar modulation. We thus ignore the flux under such low energy. Also, we have checked that the 
spectrum of 7-ray from pion decay is almost unchanged. See later discussion and Appendix. 
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the PAMELA observation at 95 % CL. We also reevaluate positron fraction and 
total (e + + e~) flux taking / y p = 2.52 above 4 GeV. The results are also shown in 
top two panels in Fig. [TJ (blue lines). One can see that the fluxes of e are almost 
unchanged. Therefore, p and p fluxes in the decaying gravitino scenario can be in a 
reasonable agreement with the PAMELA data if m 3 / 2 ~ 200 GeV. On the contrary, 
when the gravitino mass is larger, the spectrum of p/p becomes inconsistent with 
the observation of PAMELA. 

We also simulate 7-ray flux, since it is inevitably produced from cosmic-ray e ± 
via IC and bremsstrahlung processes inside Galaxy as well as via IC process in extra- 
Galactic region. Furthermore, in the present case, 7-rays are also directly produced 
as a consequence of gravitino decay and following cascade decays. In the bottom 
left panel in Fig. [TJ 7-ray from extra-Galactic region due to the decaying gravitino is 
given. Again, we take model parameters as in top two panels. Extra- Galactic 7-ray 
background (EGB) data given by EGRET [75] and Fermi-LAT [75] experiments are 
also shown. One can see that 7-ray flux from the decaying dark matter is expected to 
be smaller than the observed value when m 3 / 2 < a few TeV. On the other hand, the 
flux becomes too large to be consistent with the observation in high energy region 
Ery > 10 GeV when m 3 / 2 ~ 4 TeV; in such an energy region, the 7-ray is mainly 
from the IC process of e ± produced by the decay of dark matter. The gravitino LSP 
in this scenario with 7713/2 > 4 TeV is thus disfavored. 

In addition to extra- Galactic 7-rays, we have simulated the contribution from 
the decaying gravitino LSP inside the Galaxy. Here, we have averaged IC-induced, 
bremsstrahlung-induced, and primary 7-ray in the region b > 10° with b being 
Galactic latitude, and added them to the flux from extra-Galactic region to obtain 
the total flux from the gravitino LSP. (We have checked that the contribution from 
the bremsstrahlung process is subdominant.) The result is shown on the bottom 
right panel in Fig. [1] with the observational data of the Fermi-LAT experiment [76] . 
From the figure, one can see that the calculation leads to almost the same conclusion 
as that from the extra-Galactic 7-ray; both primary and IC-induced one are expected 
to be smaller than the observation as far as m 3 / 2 < a few TeV. Therefore, by taking 
into account the constraints from anti-proton and 7-ray, we can conclude that the 
PAMELA anomaly can be solved with the unstable gravitino LSP with m 3 / 2 ~ 
200 GeV without conflicting with other observations. 

Next, we consider the case where the gravitino LSP decays into the second gener- 
ation lepton. We have found the allowed region of 200 GeV < m 3 / 2 < 500 GeV and 
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Figure 2: Cosmic-rays from the decaying gravitino LSP in cases (1-5) — (1-8) in Table 
2J Each panel shows cosmic-ray as the same way as in Fig. [1] Results with 7 P = 2.42 
(corresponding to case (1-5), (1-6), and (1-7), where gravitino mass is taken as m.3/2 = 
250 GeV, 1 TeV, and 4 TeV, respectively) in magenta from left to right, and one with 
7 P = 2.52 (corresponding to the case (1-8) where the mass is taken to be 250 GeV) in blue 
are shown. 
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m 3/2 ^ 1 TeV. Model parameters giving the best-fit results to cosmic-ray observa- 
tions are summarized in Table H] as the cases (1-5) to (1-8). The predicted cosmic-ray 
fluxes using the best-fit parameters are shown in Fig. [2j It can be seen that the 
positron fraction as well as the total (e + + e~) flux well agree with the results of 
PAMELA and Fermi-LAT experiments irrespective of the gravitino mass (top two 
panels). It is also found that the gravitino LSP with its mass over 1 TeV results in 
a better agreement with the Fermi-LAT data than the previous case. 

As in the previous case, the observations of anti-proton flux give stringent con- 
straint on the gravitino mass. In middle two panels in Fig. HI we show the anti-proton 
flux and p/p ratio. The numerical result shows that large gravitino mass is disfa- 
vored from the observation of p/p at the PAMELA experiment. However, when the 
gravitino mass is small, the simulated p/p is consistent with the latest observation 
at 95 % C.L.. (For m 3 / 2 = 250 GeV, for example, we found Xp/ P — 31.3.) We also 
plot 7-rays in bottom two panels in Fig. [2J from which we see that the 7-ray flux 
becomes comparable to the observation when 1713/2 ~ 4 TeV. This 7-ray is mainly the 
primary one directly produced by the decay of the gravitino LSP. The IC-induced 
7-ray is not significant, because the primary e produced by the gravitino LSP is 
softer than the previous case, so that the IC-induced 7-ray is suppressed. There- 
fore, we conclude that PAMELA anomaly can be solved in the present scenario for 
200 GeV < m 3/2 < 500 GeV. 

3.3 Gravitino dark matter with tri-linear RPV 

When the gravitino LSP decays mainly through the tri-linear RPV operator given in 
Eq. (12.131) . the final state of the decay is composed of only (three) leptons. In such a 
case, constraints from anti-proton flux as well as the ratio between anti-proton and 
proton are irrelevant in the study of the decaying dark matter scenario. 

With the Xe" ariar y s i s > we have found the region consistent with e observations 
at 95% C.L. in 200 GeV < m 3/2 < 400 GeV and 1 TeV < m 3/2 < 3 TeV. (In 
the numerical calculation, we take mj = 1.2 772,3/2.) The model parameters for 
the best-fit results are again summarized in Table H] as cases (II-l) to (II-3), and 
simulated results for cosmic-ray e ± and 7-rays are shown in Fig. In upper two 
panels, positron fraction and total (e + + e~) flux are shown. Here, we consider the 
case where the component A123 is dominant compared to others. It can be seen that 
the fitting to the Fermi-LAT data becomes worse in the large 7773/2 region though it 
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Figure 3: Cosmic-rays from the decaying gravitino LSP in the cases (II-l) to (II-3) in Table 
[H where positron fraction (upper left), total (e + + e~) flux (upper light), extra-Galactic 
7-ray flux (lower left), and the flux from inside and outside of the Galaxy (lower right) 
are shown. Line type is assigned as the same as in Fig. [H i.e., signal from the decaying 
dark matter and astrophysical BG are depicted by dotted and dashed lines, respectively, 
while solid lines correspond to signal + BG. In each panel, numerical results correspond 
to (II-l), (II-2), and (II-3) from left to right (bottom to top in lower left panel), where 
gravitino mass are taken as to 3 / 2 = 250 GeV, 1 TeV, and 4 TeV, respectively. 
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is possible to explain the PAMELA data well. 

We have calculated 7-ray flux. Here, we consider the case in which A123 is the only 
relevant RPV parameter, so we also calculate primary 7-ray flux from the decay of 
r, along with the one from the IC process. The results are shown in lower two panels 
in Fig. [31 the left one is extra-Galactic contribution and the right one is extra- plus 
inner- Galactic contribution. For the extra-Galactic one, the simulated flux is much 
smaller than the observation unless 1713/2 > 4 TeV. The observation of total 7-ray 
(for b > 10°) gives almost the same upper bound on the mass. As a consequence, 
the parameter region favored by the observation of e (i.e., 200 GeV < 7723/2 J$ 400 
GeV and 1 TeV < 7723/2 < 3 TeV) are not excluded by the 7-ray observation at the 
Fermi-LAT experiment. It can be also seen that the total 7-ray flux at high energy 
region is noticeable when 7713/2 > 4 TeV. This intense flux comes from the decay of 
r. 

In addition to the case where the component Ai ?3 do minates, we have also consid- 
ered the case in which only Am or A122 is relevant P 13 l In the former case, however, 
injected e is too hard to be consistent with observations at 95% C.L. in the entire 
range of 777.3/2- On the other hand, in the latter case, the same conclusion as in the 
case A123 dominates holds except that the allowed region at the TeV scale disappears 
because the predicted total (e + + e~) flux hardly agrees with the Fermi result. 



3.4 Sneutrino dark matter with tri-linear RPV 

So far, we have discussed the gravitino LSP cases. However, there are other candi- 
dates for LSP. In particular, as we have mentioned, the sneutrino LSP may provide 
decay modes different from the gravitino LSP. In particular, with the superpotential 
given in Eq. (I2.13p . the sneutrino LSP may decay into two charged leptons. We then 
expect different cosmic-ray spectra from those of the gravitino LSP. In this sub- 
section, we briefly comment on the behaviors of cosmic-ray fluxes in the sneutrino 
dark matter scenario with tri-linear RPV. The important effects are on and 7-ray 
fluxes, but not on anti-proton flux. 

In Fig. HI we show the positron fraction (upper left) and the total (e + + e~) 
flux (upper right) for the case where the sneutrino dominantly decays through the 
process v — > fi + fi~. The parameters to give those results are summarized in Table H] 
as cases (III-l) to (III-3). 

# 13 In those cases, there is no 7-ray primarily produced by the decay. 
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Figure 4: Cosmic-rays from the decaying sneutrino LSP in the cases (III-l) to (III-3) in 
Table H3 The panel position and line contents are the same as in Fig. [3j In each panel, 
lines from left to right (bottom to top in lower left panel) correspond to the case (III-l), 
(III-2), and (III-3), where sneutrino mass is taken as rriy = 400 GeV, 1 TeV, and 4 TeV, 
respectively. 
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With the use of the Xe" ana ly s i s ! we have found 95% C.L. allowed region in 
300 GeV < mo < 500 GeV, with neglecting the possible systematic error from 
the energy calibration in the Fermi-LAT data. For niy > 500 GeV, the total flux 
becomes inconsistent with Fermi-LAT data, while positron fraction still agrees with 
the PAMELA data. If we take account of the effect of systematic error in Fermi- 
LAT data, the allowed region may become larger; in particular, adopting the 10 % 
reduction of the energy, the region 1 TeV < m^, < 3 TeV becomes allowed at 95% 
C.L. This result is consistent with [35J. 

We also give 7-r ay fl ux in Fig. HI The only relevant process to produce 7-ray here 



is the IC scattering P 14 ! In the figure, IC-induced 7-ray from extra-Galactic region 
(inner-plus extra- Galactic region) is given lower left (right) panel. Although the 
flux is less than the observed data in b > 10°, IC-induced 7-ray in the extra-Galactic 
region becomes comparable to or larger than the data when > 4 TeV. This is the 
same result obtained in the previous three-body decay case. We thus conclude that 
the parameter region 300 GeV < < 500 GeV gives good fit with PAMELA and 
Fermi-LAT (e + + e~) data and the region is not constrained by 7-ray observation. 

We have also considered the decaying scenario of v — > e + e~. In such a case, 
however, we could not find the allowed region. This is because the monochromatic 
electron and positron gives very sharp edge in the flux so that the flux does not agree 
with the data of (e + + e~) flux though the fit with the PAMELA data is good. 



4 Conclusion 

In this article, we have calculated cosmic-ray fluxes from decaying dark matter as 
well as background in the same propagation model by using GALPROP. Aiming 
for explaining the PAMELA anomaly with being consistent with other cosmic-ray 
observations, we have reevaluated the background cosmic-ray fluxes. Cosmic-rays 
from decaying dark matter, on the other hand, is strongly dependent on the spectra 
of final-state particles. If one specifies the distributions of final-state particles, the 
cosmic-ray spectra are determined independently of detailed framework of the model 
of decaying dark matter. 

To make our discussion concrete, we have studied gravitino dark matter in R- 

# 14 There could exist contribution from final state radiation (FSR). FSR may give hard 7-ray and 
also give important contribution when the leptonic decay is suppressed by kinematics or chirality; 
however it is not such a case which we consider. 
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parity violated supersymmetric model. Under i?-parity violation, gravitino domi- 
nantly decays to Wli m bi-linear RPV, while it decays to Vil^lJ in tri-linear one. 
In the former scenario, we have found that the simulated cosmic-ray e ± flux from 
dark matter and background agrees with the PAMELA and Fermi-LAT data, ir- 
respective of gravitino mass. However, it has been shown that the production of 
cosmic-ray p becomes enhanced when the mass is large; thus the gravitino mass 
larger than ~ 300 GeV is disfavored. For 7-ray, the flux is consistent with the ob- 
servation as far as m 3 / 2 < 4 TeV. In the latter case, it has been found that the 
PAMELA anomaly can be explained in the mass region 200 GeV < m 3 / 2 < 400 GeV 
or 1 TeV < m 3 / 2 < 3 TeV, being consistent with Fermi-LAT data. When the mass is 
larger than ~ 4 TeV, the fit with each data becomes worse. In addition, IC-induced 
7-ray constrain such large mass region. 

We have also considered the case where sneutrino is the LSP, assuming that it 
decays into lepton pair via tri-linear RPV interaction. We have seen that, when the 
dominant decay mode is v — > fJ> + fi~, the positron fraction can be in a good agree- 
ment with the PAMELA data without conflicting the Fermi result when 300 GeV < 
rriy < 500 GeV. With such a choice of the sneutrino mass, 7-ray flux induced by the 
dark matter decay is much smaller than the observed one. In addition, in this case, 
it should be noted that the constraint from anti-proton flux is irrelevant because 
hadrons are hardly produced by the decay of v. 

Acknowledgments 

This work was supported in part by Research Fellowships of the Japan Society for the 
Promotion of Science for Young Scientists (K.I.), and by the Grant-in- Aid for Scien- 
tific Research from the Ministry of Education, Science, Sports, and Culture of Japan, 
No. 21740174 (S.M.) and Nos. 22540263 and 22244021 (T.M.). We also thank orga- 
nizers, especially, at Niigata University of SI2010(YITP-W-10-07) in Fuji-Yoshida, 
and the Yukawa Institute for Theoretical Physics at Kyoto University. Discussions 
during the workshop were useful to finalize this project (S.M. and T.M.). 



21 



Background proton 



1 1 

BESS > 
CAPRICE 
ATIC2 i 



T (GeV) 




1000 



0.4 
0.35 

0.3 
0.25 

O 

m 0.2 
0.15 | 

0.1 
0.05 




0.01 



Boron to Carbon ratio 



T 



HEAO-31- 

ATIC2 
CREAM h 

I . . 



H 6 



0.1 



1 10 

T (GeV) 



100 



1000 



Figure 5: Background proton flux (the left panel) and the B/C ratio (the right panel). 
Here, the blue (magenta) line shows the simulated result by 7 P = 2.52 (2.42) above 4 GeV. 
Other parameters are taken as the same in Conventional model. 



A Background Cosmic-Rays 

We consider background cosmic-rays predicted in the Conventional model, except for 
taking 7 P = 2.52 (instead of 2.42) above 4 GeV. The main purpose in this appendix 
is to show that the change of 7 P does not significantly affect background cosmic-ray 
fluxes, and that such a choice is consistent with observations. Cosmic-ray p flux and 
the B/C ratio are shown in Fig. 5, where simulated results are depicted with a blue 
line. The results with 7 P = 2.42 above 4 GeV are also shown for comparison with 
a magenta line. For the p flux shown in the left panel of the figure, it can be seen 
that the spectrum with 7 P = 2.52 becomes slightly softer than that wit h 7 P = 2.42 



#15 



as expected, and both spectra are well consistent with observations. 

On the other hand, the B/C ratio is expected to be hardly changed. This is 
because the fluxes of B and C are mostly determined by the secondary and primary 
fluxes from SNR and hence the B/C ratio depends strongly on diffusion parameters 
but is insensitive to SNR spectra. In fact, as shown in the right panel of Fig. \K\ 
there is no difference between 7 P = 2.42 and 2.52 cases, which are consistent with 
observations. (The intensity in E < 1 GeV is larger than the one appeared in 
Ref. [32] . This is due to the different choice of solar modulation potential; we take 
= 550 MV, whereas = 450 MV is chosen in Ref. [32].) 



# 15 In the region E < 10 GeV, the flux with 7 P = 2.52 seems to slightly exceed the observations, it 
can be, however, optimized by choosing a proper normalization of the SNR p flux. Moreover, in this 
region, the flux is affected by the effect of solar modulation, which leads to further uncertainties. 



22 



References 

[1] O. Adriani et al. [PAMELA Collaboration], Nature 458, 607 (2009). 

[2] D. Hooper, P. Blasi and P. D. Serpico, JCAP 0901, 025 (2009). 

[3] Y. Nomura and J. Thaler, Phys. Rev. D 79, 075008 (2009). 

[4] P. f. Yin, Q. Yuan, J. Liu, J. Zhang, X. j. Bi and S. h. Zhu, Phys. Rev. D 79, 
023512 (2009). 

[5] Y. Bai and Z. Han, Phys. Rev. D 79, 095023 (2009). 

[6] C. R. Chen, F. Takahashi and T. T. Yanagida, Phys. Lett. B 673, 255 (2009). 

[7] K. Hamaguchi, E. Nakamura, S. Shirai and T. T. Yanagida, Phys. Lett. B 674, 
299 (2009). 

[8] E. Ponton and L. Randall, JHEP 0904, 080 (2009). 

[9] A. Ibarra and D. Tran, JCAP 0902, 021 (2009). 

[10] C. R. Chen, M. M. Nojiri, F. Takahashi and T. T. Yanagida, Prog. Theor. Phys. 
122, 553 (2009). 

[11] A. Arvanitaki, S. Dimopoulos, S. Dubovsky, P. W. Graham, R. Harnik and 
S. Rajendran, Phys. Rev. D 79, 105022 (2009). 

[12] K. Hamaguchi, S. Shirai and T. T. Yanagida, Phys. Lett. B 673, 247 (2009). 

[13] I. Gogoladze, R. Khalid, Q. Shan and H. Yuksel, Phys. Rev. D 79, 055019 
(2009). 

[14] K. Hamaguchi, F. Takahashi and T. T. Yanagida, Phys. Lett. B 677, 59 (2009). 
[15] E. Nardi, F. Sannino and A. Strumia, JCAP 0901, 043 (2009). 
[16] C. R. Chen and F. Takahashi, JCAP 0902, 004 (2009). 

[17] D. Aristizabal Sierra, D. Restrepo and O. Zapata, Phys. Rev. D 80, 055010 

(2009) . 

[18] D. A. Demir, L. L. Everett, M. Frank, L. Selbuz and I. Turan, Phys. Rev. D 
81, 035019 (2010). 

[19] L. Zhang, C. Weniger, L. Maccione, J. Redondo and G. Sigl, JCAP 1006, 027 

(2010) . 

[20] C. D. Carone, J. Erlich and R. Primulando. larXiv: 1008.06421 [hep-ph]. 



23 



[21] K. Ishiwata, S. Matsumoto and T. Moroi, Phys. Lett. B 675, 446 (2009). 

[22] M. Cirelli, M. Kadastik, M. Raidal and A. Strumia, Nucl. Phys. B 813, 1 (2009). 

[23] I. Cholis, D. P. Finkbeiner, L. Goodenough and N. Weiner, JCAP 0912, 007 
(2009). 

[24] D. Feldman, Z. Liu and P. Nath, Phys. Rev. D 79, 063509 (2009). 

[25] P. J. Fox and E. Poppitz, Phys. Rev. D 79, 083528 (2009). 

[26] L. Bergstrom, T. Bringmann and J. Edsjo, Phys. Rev. D 78, 103520 (2008). 

[27] V. Barger, W. Y. Keung, D. Marfatia and G. Shaughnessy, Phys. Lett. B 672, 
141 (2009). 

[28] A. E. Nelson and C. Spitzer. lirXiv:0810.5167 [hep-ph]. 

[29] R. Harnik and G. D. Kribs, Phys. Rev. D 79, 095007 (2009). 

[30] M. A. DuVernois et al, Astrophys. J. 559 296 (2001). 

[31] K. Ishiwata, S. Matsumoto and T. Moroi, Phys. Rev. D 78, 063505 (2008). 

[32] A. Ibarra and D. Tran, JCAP 0807, 002 (2008). 

[33] GALPROP Homepage, |http :"/ /galprop . Stanford . edu/ 



[34] J. Zhang, X. J. Bi, J. Liu, S. M. Liu, P. F. Yin, Q. Yuan and S. H. Zhu, Phys. 
Rev. D 80, 023007 (2009). 

[35] V. Barger, Y. Gao, W. Y. Keung, D. Marfatia and G. Shaughnessy, Phys. Lett. 
B 678, 283 (2009). 

[36] V. Barger, Y. Gao, W. Y. Keung and D. Marfatia, Phys. Rev. D 80, 063537 

(2009) . 

[37] J. H. Huh and J. E. Kim, Phys. Rev. D 80, 075012 (2009). 

[38] N. E. Bomark, S. Lola, P. Osland and A. R. Raklev, Phys. Lett. B 686, 152 

(2010) . 

[39] T. Lin, D. P. Finkbeiner and G. Dobler. larXiv: 1004.09891 [astro-ph.CO]. 

[40] R. C. Cotta, J. A. Conley, J. S. Gainer, J. L. Hewett and T. G. Rizzo, 
larXiv:1007.5520l [hep-ph]. 

[41] A. W. Strong, I. V. Moskalenko and O. Reimer, Astrophys. J. 537, 763 (2000) 
[Erratum-ibid. 541, 1109 (2000)]. 

[42] A. W. Strong, I. V. Moskalenko and O. Reimer, Astrophys. J. 613, 962 (2004). 

24 



[43] A. A. Abdo et al. [The Fermi-LAT Collaboration], Phys. Rev. Lett. 102, 181101 
(2009). 



[44; 
[45 

[46 

[47 

[48 
[49 
[50 
[51 
[52 
[53 
[54 
[55 
[56 
[57; 



L. J. Gleeson and W. I. Axford, ApJ 154, 1101 (1968). 

T. A. Ensslin, P. L. Biermann, P. P. Kronberg and X. P. Wu, Astrophys. J. 
477, 560 (1997); 

A. Loeb and E. Waxman, Nature 405, 156 (2000); F. Miniati, Mon. Not. Roy. 
Astron. Soc. 337, 199 (2002). 

Y. T. Gao, F. W. Stecker, M. Gleiser and D. B. Cline, Astrophys. J. 361, L37 
(1990); A. Dolgov and J. Silk, Phys. Rev. D 47, 4244 (1993). 

J. F. Navarro, C. S. Frenk and S. D. M. White, Astrophys. J. 490, 493 (1997). 

K. Ishiwata, S. Matsumoto and T. Moroi, Phys. Lett. B 679, 1 (2009). 

M. Cirelli, P. Panci and P. D. Serpico, Nucl. Phys. B 840, 284 (2010). 

G. Hutsi, A. Hektor and M. Raidal, JCAP 1007, 008 (2010). 

J. M. Cline, A. C. Vincent and W. Xue, Phys. Rev. D 81, 083512 (2010). 

A. A. Abdo et al, Astrophys. J. 712, 147 (2010). 

M. Perelstein and B. Shakya, |arXivil002,4588 [astro-ph.HE]. 



M. R. Buckley and D. Hooper, |axXiv: 1004. 1644 [hep-ph]. 

R. Essig, N. Sehgal and L. E. Strigari, Phys. Rev. D 80, 023506 (2009). 

R. Essig, N. Sehgal, L. E. Strigari, M. Geha and J. D. Simon, larXiv:1007.4T99l 
[astro-ph.CO]. 

[58] M. Perelstein and B. Shakya, arXiv:1007.00T8l [astro-ph.HE], 

[59] T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605, 026 (2006). 

[60] T. Moroi, H. Murayama and M. Yamaguchi, Phys. Lett. B 303, 289 (1993). 

[61] J. L. Feng, A. Rajaraman and F. Takayama, Phys. Rev. Lett. 91, 011302 (2003). 

[62] K. Ishiwata, S. Matsumoto and T. Moroi, Phys. Rev. D 77, 035004 (2008). 

[63] T. Asaka, K. Ishiwata and T. Moroi, Phys. Rev. D 73, 051301 (2006). 

[64] B. A. Campbell, S. Davidson, J. R. Ellis and K. A. Olive, Phys. Lett. B 256, 
484 (1991); W. Fischler, G. F. Giudice, R. G. Leigh and S. Paban, Phys. Lett. B 
258, 45 (1991); H. K. Dreiner and G. G. Ross, Nucl. Phys. B 410, 188 (1993). 

[65] L. J. Hall, T. Moroi and H. Murayama, Phys. Lett. B 424, 305 (1998). 

25 



[66] K. Ishiwata, S. Matsumoto and T. Moroi, JHEP 0905, 110 (2009). 

[67] J. Chang et ai, Nature 456, 362 (2008). 

[68] S. Torii et ai, arXiv: 0809 .07601 [astro-ph], 

[69] D. Hooper and K. M. Zurek, Phys. Lett. B 691, 18 (2010). 

[70] S. Orito et al. [BESS Collaboration], Phys. Rev. Lett. 84, 1078 (2000); 
Y. Asaoka et al, Phys. Rev. Lett. 88, 051101 (2002); K. Abe et al, Phys. 
Lett. B 670, 103 (2008). 

[71] M. Boezio et al. [WiZard/CAPRICE Collaboration], Astrophys. J. 561, 787 
(2001). 

[72] PAMELA Collaboration. larXiv: 1007.08211 [astro-ph.HE]. 

[73] F. Donato, N. Fornengo, D. Maurin and P. Salati, Phys. Rev. D 69, 063501 
(2004). 

[74] O. Adriani et ai, Phys. Rev. Lett. 102, 051101 (2009). 

[75] P. Sreekumar et al. [EGRET Collaboration], Astrophys. J. 494, 523 (1998). 

[76] A. A. Abdo et al. [The Fermi-LAT collaboration], Phys. Rev. Lett. 104, 101101 
(2010). 



26 



